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ABSTRACT 


In this thesis, the algorithms for an A-7E aircraft perfor- 
mance calculator were developed and then implemented on three 
small data processors of different programming levels and storage 
capabilities. 

The utility of data is a function of several variables in- 
cluding accuracy and availability. The problem of retrieving 
me@eermeance data from the Naval Air Training and Operating Pro- 
cedures Standardization (NATOPS) Manuals is significantly les- 
sened by the devices demonstrated in this investigation. Nine 
performance chart groups, ylelding data usually considered 
necessary for flight, were reduced to a series of analytical 
expressions. These analytical expressions were demonstrated to 
reproduce NATOPS Manual data to a high degree of accuracy. 

Implementation was demonstrated on a desk computer, a hand 


held calculator and a microprocessor. 
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i SE ObUCT LON 


The Naval Air Training and Operating Standardization 
(NATOPS) Manual is the official standard of the United States 
Navy for "...information on all aircraft systems, performance 
data, and operating procedures required for safe and effective 
operations." [1] 

The purpose of this thesis was to develop algorithms of 
the more often used NATOPS performance charts for the A-7E 
alrcraft, examine their accuracy and implement them on small 
data processors that might be adaptable to shipboard or aircraft 
onboard use. The interpretation of NATOPS performance charts 
1s an error prone and time consuming procedure even for experi- 
enced users. The need for a system to eliminate this laborious 
process has been fully documented in a thesis completed in June 
mayo by LCDR W.M. Siegel [2]. In his investigation, LCDR Siegel 
Sevised an efficient procedure to develop algorithms from the 
NATOPS performance charts and exercised this procedure on the 
problems of "Takeoff Ground Roll Distance" and "Takeoff Airspeed". 

This investigation is an extension of the aforementioned 
work. The original scope of this investigation was to develop 
algorithms for eleven of the most often used performance problem 
chart groups and implement them on the Texas Instruments-59 
(TI-59) hand held calculator (HHC). All of the NATOPS perfor- 
mance charts were not reduced because of research time 


limitations. Of the eleven performance chart groups studied, 





two performance problems, "Time to Climb" and "Fuel Required 

to Climb" were rejected because of implementation difficulties 
on the TI-59 HHC (discussed fully in "Development Difficulties"). 
Therefore, nine performance chart groups were reduced to analy- 
tical expressions and implemented on the TI-593 HHC. To show 
further possibilities and feasibility of implementation of the 
algorithms, they were 1) fully implemented on the Hewlett 
Packard-9830 (HP-9830) desk computer, 2) demonstrated on a 
microprocessor (INTEL Corporation Microcomputer System-48), 

and 3) considered for implementation on the A-7E onboard digi- 
tal computer and a microprocessor utilizing a recently developed 
number processing chip by the National Semiconductor Corporation 


m@ai57109). 


IL{a 





ie Sey CORMENT 


fee GUIDELINES 

The scope of this investigation was established after a. 
firm set of guidelines was defined. 

Being the official United States Navy standard for the 
A-7E aircraft, the A-7E NATOPS Manual was the sole source of 
performance data used to develop the algorithms. As such, 
amd being subject to changes during the aircraft's life cycle, 
the need for possible future updates to the algorithms was 
acknowledged. The effective date of the NATOPS Manual from 
which these algorithms were developed 1s March 1975. Since the 
performance data yielded by the algorithms was identical to 
NATOPS Manual performance curves, the same restrictions and 
limitations apply. For example, takeoff airspeed calculation 
restricts the NATOPS Manual user to trailing edge flap positions 
between 20 and 40 degrees down (Figure 1). For that reason, one 
could not expect to calculate the flaps up takeoff airspeed 
using the developed algorithms. An additional feature provided 
by the algorithms was higher order interpolation. While the 
Benexperienced NATOPS user might attempt to interpolate linearly 
between non-linearly oe curves, the algorithms do not. 

An important guideline for the user's benefit was to ensure 
the execution of these algorithms after implementation was 
Simple enough so very little training was required for the 


users. Intended users were Naval Flight Officers and Aviators. 


Jae 





ENGINE: TF41-A-2 
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MILITARY RATEO THRUST 
LANOING CONFIGURATION 
LEAOING EDGE FLAPS DOWN 


CONOITIONS: 


TAKEOFF SPEED (A-7E) 


MOOEL: A-7E 
DATE. 


NAVAIR 01-45AAE-1 


44 
76E€ 254-02- 72 










36 












Takeoff Speed NATOPS Chart 














# ~~ = 
~ © row) 
aS 
7 ; 
x6 ne 
X35 = 
HN ¥ 
Qe 
2 3 2 
2 3 2 = 
a Oo 
Sea Y 
ae 
J 3 < 
2 3 00 
ne ; a 
“ i 3 
aah oF. 
Aah (6) as 
eas e2 
ee 2 a 
3 5 ce) oe a 
pate 2: | (Bet eae eae 
-— = ie 
LS eC aE searaNe 
JSC CC SeNSE Ge 
™N 
= : a 
e e i] w 
Vee eee eae ee 
rc 
\ O 
\ wd 
SVS a ee 
ib ° 1 
” 
Reese eM A TA 2 
\ \ . : 
° uw 
> e 
uw 
2 o = 8 re 
~ ~- ~ R = 


SVIN — ALIDOTZA 


11-2) 


12 





Not included in the scope of this thesis is an introduc- 
tion to the TI-59 HHC, HP~9830 desk computer and the INTEL 
Microcomputer System-48; however, to follow the computer pro- 
grams written for these devices would required their basic 
understanding. 

Another guideline established was that the performance 
calculators be light and small enough to be physically suited 
for its environment. For example, the TI-59 calculator and 
microprocessor could be used in a cockpit, briefing room or 
Air Operations Center. The HP~9830 desk computer would be 
restricted from cockpit use. 

Reliability was a necessary guideline. 

To make algorithm implementation on the TI-53 HHC feasible 
amd since the program storing chip, the Continuous Read Only 
Memory (CROM), was limited to 5000 calculator program steps, 
the library of nine programs was required to fit into that 
space [3]. 

Finally, accuracy was a necessary consideration. The 
results obtained from the algorithms were required to be at 
least as accurate as following the performance charts manually. 
These accuracy requirements established were: One knot of air- 
Meee 100 feet of altitude or ground roll distance, 100 pounds 
of weight, ten seconds of time and one nautical mile cf dis- 


Tance. 
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Pee PERFORMANCE CHART REDUCTION 

The reduction of the NATOPS Manual performance curves into 
analytical expressions was accomplished by a historically proven 
mathematical procedure, "least squares curve fitting". This 
method was applied to certain A-7E performance data by LCDR 
W.M. Siegel (see Introduction, Section I). His brief explana- 
tion of the "Least Squares Fit Approximation (LSFA)" is included 
in Appendix A. 

Many performance charts from the NATOPS Manual contain 
three variables (two independent, one dependent) and are de- 
picted as a two-dimensional space with the third dimension 
illustrated by a family of curves. The reduction of such a 
chart can be accomplished as follows: 

1. Determine order of curves in family (1.e, second order, 
= A, + Aox + A,x*). 

2. Apply LSFA to every member of the family of curves. 

3. Since the order of the curve families may vary, a 


general curve family could be depicted as follows: 


= Zz 
y = Ay si Ay 5X i: A, 3% are ine Ay me (for curve Z4) 
a 2 n-l 
y = Any + An o® + An3% fears AynX (for curve Zn) 
: = i : ; 2 , n-L : : 
y = Anl + Ano% 1 A_43% : ae ae x (for curve a) 


pee apply LSFA to the coefficients. For example, plot Ai. 


Any > 5 oe Anal versus 2 


s 2 
A. = Bay + B52 + B42 t -++ Bs, 


Doing the same with all coefficients, 


1? Zo? respectively, yielding 


r=-1 
Z: J 


14 





2 r-1 


Hy S20, ay aie ie Soe 
oa : ; ; ; on 
AL = Bal F B22 # B32 + -++-B 2 
5. Given z and x, y can now be calculated by: 
a. Computing coefficients from equations generated in 


Step 4. 
b. Applying coefficients to y = A, + A,X a we 
Seite 1S important to note that although all curve family 
members must be of identical order, the equations representing 
Ene coefficients as a function of "z" need not be of similar 
perder. 

Although applying LSFA to the family of curves and then 
Bemtheir coefficients was the normal method of chart reduction, 
it was not always used for the following reasons: 

a. Some charts were two-dimensional (LSFA still used). 

b. Some charts were reduced by inspection. 

(1) Linear curve families with linear spacing. 
(2) Time, distance, speed charts (d = v/t). 

c. Algorithm anomalies (see "Development Difficulties"). 

When used, the LSFA was accomplished by a program pre- 
written by the Hewlett Packard Corporation for use with the 
Me=-Jo30. This program, although greatly facilitating the 
development portion of this investigation, was written for a 
two-dimensional problem and had to be executed at least once 
for each curve and once for each set of coefficients. 


A disting of all of the equations making up the perfor- 


Mance algorithms are contained in Appendix C. The A-7E 
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performance chart groups from which they were developed are 
contained in Appendix B. They are in order: 

1. Low Level Cruise Performance. 

2. Takeoff Ground Roll Distance. 

3. Maximum Range Cruise Time and Speed at Constant 
Sepeatude. 

4. Maximum Range Cruise Fuel Required at Constant 
Altitude. 

5. Maximum Range Climb Airspeed Schedule. 

6. Takeoff Airspeed. 

7. Maximum Refusal Airspeed. 

ca Cptimumeendurance Altitude. 

9. Cruise Ceiling. 

Future reference in this thesis is made to algorithms and 


programs by the numbers above. 


See EXAMPLE OF CHART REDUCTION 

An example of the procedure discussed in the previous sec- 
tion is presented below. The chart chosen for reduction is the 
lower graph of Figure 2, from Phase II of the A-7E Cruise Per- 
formance chart group. 

Eyeinspection, all Ay and A. coefficients are equal to 
zero. The curves appear parabolic and therefore second order, 
yielding N = AjM*. The example follows: 


N imeawmnecd late Pes i 


M mach number 


D = drag count 
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CRUISE PERFORMANCE (A-/E} 


PHASE II — AIRCRAFT REFERENCE NUMBER 


MODEL: A-7E ENGINE: TF41-A-2 
DATA BASIS: FLIGHT TEST FUEL GRADE: JP-5 
DATE: NOVEMBER 1971 FUEL DENSITY: 6.8 LB/GAL 
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DRAG COUNT LINE CURVE EQUATION 


50 NS UU. Se EB IGE 
100 Mea 2. Tae 7M- 
150 N = 4.1658M? 
200 N = 5.5530M? 
300 Neg=ees 327 aM 
u00 = iene. 2M = 


Byeplotting the A, coefficients versus D (drag count), the 


3 
BerA yields: 
A, = $Gteo7 sZe=3) + .027743DRand therefore, 
N = G(4.3732E-3) + (.027743D)M*. 
This was a particularly simple ataauare to reduce but illus- 


trates the procedure. 


Meee DEVELOPMENT DIFFICULTIES 

The normal method of reducing performance curves did net 
always yleld useful information. One reason was although the 
NATOPS Manual Performance curves were constructed from experi- 
mental data, families of curves occasionally had very unusual 
Spacing. They also were not always a true curve family; that 
is, they were of varying order. This can be visually detected 
in the lower graph of Phase III of the A-7E Cruise Performance 
Chart group (Figure 3). The unequal and varying spacing be- 
tween curves with different "reference numbers" is obvious. 
Although the coefficients for each curve can be calculated, 
the coefficients determined for a LSFA equation for an inter- 


memes Curve would-be incorrect. To be usable for the normal 


lis 
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method of chart reduction, a chart must have equal, constantly 
increasing, or constantly decreasing spacing between curves. 
When such an incompatible chart was encountered, it was neces- 
sary to interpolate between them. Two chart groups eliminated 
from consideration, "Fuel Required" and "Time to Climb from Sea 
Level to Selected Altitude", contained so many such curves 
(11), that very high order expressions would have been re- 
quired to compute the coefficients, making implementation on 
the TI-59 HHC impractical. The A-7E Cruise Performance lower 
chart of Phase I had the same anomaly (Figure 4). Because of 
the importance of the low level mission, however, the algorithm 
for this chart was developed, for sea level only though. The 
multiple algorithm was not developed but could have been for 
implementation on a desk computer. 

Another reason a straight application of LSFA was not 
always appropriate was the uniqueness of the upper graph of 
Phase I of the A-7E Cruise Performance chart group (Figure 4). 
This chart requires entry from the lower chart. A line is 
traced upward until the user contacts the appropriate Drag 
Count Line (dotted lines). The first pass through the Mach 
Number axis, a result of the lower chart, was defined M*. 
Instead of now tracing horizontally to the Transfer Scale axis 
Menas value defined TS*), one must trace "between the solid 
guidelines" to the interception with a line traced vertically 
upward from the desired Mach number, M. The Transfer Scale 


would now be manually obtained by tracing horizontally to the 
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Cruise Performance 
Phase I NATOPS Chart 


cu 








vertical axis. To develop the algorithm for this problen, 

the equations of the guidelines were also calculated as a func- 
tion of Mach number. The values of the Transfer Scale resulting 
from M* intercepting the guidelines and tracing horizontally 
forthe vertical axis were called TS +; TS, ...TS*, from top 
Membottom. The original position, (M*, TS*), could now be 


determined in relation to (M*, Ts *) and (M*, TS re oe ae and 


ao 
"nt+1" indicate the upper and lower guidelines, respectively, 
which bracket (M*, TS*). This ratio provided the initial posi- 
tion relative to the guidelines: 

=) 


R = (TS*-TS ®)/CTS8-TS 
iar Jb n 


+ ila JL 


Using the desired Mach number, M, the Transfer Scales for 
the same two enclosing guidelines were calculated (TS. and 
TS o41?° The final position relative to the guidelines was 


maintained using the original ratio by solving: 


m= (x-TS VAG Ss= LS rosa on 
mr n n 


+ al 


"x" as the Transfer Scale with which the user now proceeds to 
Phase III of this performance chart group. Figure 5 depicts 


this problem graphically. 


me ACCURACY 

A large number of results comparisons between the generated 
algorithms and manually traced performance problems were made. 
An infinite number of comparisons would be required to check 
all possibilities, but since the mathematical theory was so 
basic, the number of checks accomplished were considered suf- 


mietent. 


Ze 





MACH NUMBER 


Figure 5 


Guideline Chart Solution 


TRANSFER SCALE 





All nine algorithms were checked for accuracy on the 
HP-9830 desk computer. The number of checks for each algo- 
rithm was proportional to the ease of manually tracing through 
the performance charts. The author spent considerable time 
obtaining performance results from the NATOPS Manual charts 
and a relatively small amount of time computing the problems 
on the desk computer once the algorithms had been implemented. 
In a significant number of instances, the results disagreed, 
but after rechecking, the solution obtained manually was in 
error. This supported the contention that manual manipulation 
of the performance charts 1S an error prone procedure, even 
with an experienced user. 

In a few rare instances, the author entered the required 
given data incorrectly into the desk computer. These miskeying 
errors, not procedural, were noticed as soon as the answer was 
produced. A user familiar with the A-7E performance character- 
istics would normally notice an answer resulting from grossly 
incorrect data input. It is acknowledged, however, that there 
is no failsafe check on the programs. When using a desk com- 
puter, the required input data can be printed along with the 
answer to ensure the user of the correctness of the input data. 
For a hand held calculator, however, computing a performance 
problem twice would provide a check, which is what many NATOPS 
Manual users often do. As with all computer programs, a desired 


result requires accurate input data. 
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Except for those noted below, the results of programs 
checked (using five significant figures) were indistinguish- 
able from the answers obtained by manually manipulating the 
performance charts. Answers produced from the algorithms 


were rounded off to the nearest digit. 


PROGRAM MAXIMUM DEVIATION 
Maximum Refusal Speed 2 knots 
Takeoff Airspeed es agtehe 


Zo 





ee ee eect A ch ON 


fee DeEoK COMPUTER 

The use of a desk computer capable of producing A-7E ee 
formance information within seconds (less than three seconds 
computation time for the longest algorithm) would be ideal for 
a squadron briefing room or Air Operations Center use. The 
HP~9830 desk computer was used for this implementation stage. 
Very little training would be required for personnel to load 
the programs stored on a cassette tape cartridge and execute 
them. 

A knowledge of "basic" computer language is required to 
fully understand the nine HP-9830 programs in Appendix D [4]. 
The nine programs are in the same order as the algorithms of 
Appendix C. 

Only in the Low Level Cruise Performance program are sub- 
routines required for linear interpolation or for the iterative 
M-enod to find the Transfer Scale (see "Development Difficulties". 
All other programs are straight forward, sequential computations. 
In these programs, the coefficients defining a curve (y = f(x)) 
for a given set of conditions are calculated. That chart 
result, "y", 1s then calculated for the given independent vari- 
able "x". The next chart of that group 1S similarly treated 
mae sO On until the "final result" is achieved. 

The HP-9830 programs are very useful since they prompt the 


user to supply the correct information. Most of the programs 
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"nequest, then accept" those inputs required for the applicable 
NATOPS Manual performance chart. The HP-9$830 then prints the 
data just entered (ensuring the user that data input was as 
desired) followed quickly by the solution. The computer is 
instantly ready to receive new data for another calculation. 

Programs 1, 2, 3, 4 and 7 (as identified in "Performance 
Chart Reduction, section II-B), are written in this "request, 
hen accept" format. The shorter programs, 5, 65, 8 and 9, 
were written with an initial set of input data already in the 
program. This format allowed the computer to step incrementally 
through the allowable range of values for the input data, thus 
calculating a "table of performance data" for the applicable 
performance chart group. These programs are easily altered to 
the "request, then accept" format by some simple edit commands 
C4]. 

The variables used in the programs are defined following 


Sa@m program in Appendix D. 


PeeeaAND HELD CALCULATOR 

The many favorable features of the hand held calculator 
encouraged its implementation of the performance algorithms. 
Its small size allowed consideration for use in the cockpit. 
Its simplicity and reliability was an advantage making it 
especially suited for users of varying experience (including 
no experience). Although its execution speed was the slowest 
of all devices used, the computation time was still much faster 


than using the NATOPS Manual. 


oa 


The Texas Instruments-59 (TI-59) programmable hand held 
calculator (HHC) was selected for implementation. This selec- 
tion was made for several reasons. At the time, it was the 
only calculator available to the author which allowed permanent 
program storage (on magnetic cards). Additionally, the Texas 
Instruments Corporation had the capability to combine all pre- 
written performance programs, up to a 5000 program step limit, 
onto a Continuous Read Only Memory (CROM) chip, making the 
A-7E performance programs a permanent part of the calculator. 
mas CROM chip can also be used on the less expensive TI-58 HHC. 
These features made the TI-58/59 (with CROM) a practical system 
form the A-7E Naval Aviation community. 

One might consider the calculator's inability to prompt the 
user for inputs a shortcoming of this implementation candidate, 
but a company spokesman, Mr. Richard Cuthbert, stated a new 
Bteemeould be fitted onto the calculator, identifying different 
buttons with the input data categories such as GW for gross 
weight, FLPS for flap position, T(°C) for temperature, and so 
on ae 

Some time was required for the author familiarization with 
the TI-59 HHC and its capabilities. For a detailed explanation 
of comments in this section involving TI-59 programming and 
Smpendix £, consult the user's manual [5]. 

All programs were entered with the calculator memory parti- 
tioned to allow 879 program steps and ten memory storage 
locations. The loss of program steps in order to provide 


SGerticient storage locations (ten to one) was the reason for 
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Meee tioning in this manner. Only five significant figures 
were considered necessary for computational accuracy. Consider- 
ing the number possibilities (1.2345 to 1.2345E-12) might take 
from six to ten program steps, this was less than the absolute 
ten program steps sacrificed for a storage location. The ten 
memory storage locations were used to store the input data at 
program execution start but were often reused after the input 
data storage was no longer required. 

The programming language level of the TI-59 HHC is below 
the HP-9830's and above a microprocessor's (discussed later) 
in sophistication. The algorithms were computed in a more 
Space-sSsaving manner than on the HP-9830. For example, in com- 
Sijeums 2 first order polynomial, the HP-9830 program functioned 
as follows: 


B(O) 


tl 


Ay, * Ay 22 


e( 1) = Any 7° Ano2 
yee BCG) + BCL). 
The TI-59 HHC was programmed to compute as follows: 
ae 7 Y- 

In the Low Level Cruise Performance program, the linear 
interpolation and iterative methods to follow guidelines (dis- 
cussed in previous section) was still accomplished using the 
more tedious TI-59 HHC language. 

Using the partitioning already described, a program limit 


of 879 program steps was imposed (filling two magnetic cards). 


Two programs, "Takeoff Ground Roll Distance" and "Low Level 
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Cruise Performance", exceeded this limit and had to be continued 
on extra cards. These programs were written to allow storage 

of an intermediate result into the T-register. The rest of the 
cards could then be read in, any lost or newly acquired input 
data entered, and program execution would continue, automatically 
retrieving the stored intermediate result from the T-register. 
These artificial necessities for program completion uSing the 
magnetic cards would not be necessary if the programs were 

stored permanently in the CROM. 

The total number of steps required for the nine performance 
algorithms programmed on the TI-59 HHC was 5461 steps. By sub- 
routining (340 steps of programming are common to two programs), 
the total number could be reduced to 5121 steps. The elimina- 
tion of the artificial steps required for the oversized programs 
would reduce the overage more. The sole intent of this imple- 
mentation phase was not to fit these nine programs into the 
2000 step CROM. If the inclusion of all nine programs was 
desired, streamlining aid offered by engineers from the Texas 
Mieeruments Corporation plus the reduction of significant 
figures in a non-critical area would accomplish this. 

The program listings, storage location usage, user instruc- 


tions, and execution times are included in Appendix E. 


mee LCROPROCESSOR 


ieeoingle Board Computer using Software for Mathematical 
Operations 


The single board computer (SBC) implementation was 


investigated both as an extension of thesis work and to meet 
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the course objectives of AE-4900, Air Data Systems. Work 
toward this effort was also done by LCDR W.M. Siegel. The 
performance algorithms were to be processed on a SBC using an 
INTEL Corporation 8048 Programmable Read Only Memory (PROM), 
external random access memory (RAM) and a program counter. 
software development was completed on the INTEL Prompt-48 
(Microcomputer System-48 language) using an INTEL 8035 arith- 
Metre Jogic unit (ALU). Although a SBC using the 8048 PROM 
and requiring a digital keyboard and display was never actually 
constructed because of the time limitations, the software opera- 
tion was successfully demonstrated on the Prompt—48. 

To preserve the programs between operation periods, 
the Prompt 48 was hand wired as specified in the user's manual 
to an ASR-35 Teletype set which allowed paper tape storage [6]. 
The Prompt-48 provided 1024 by two bytes of RAM and 64 by two 
bytes of resident memory. Although the MCS-48 instruction set 
will not be discussed in this thesis, a basic understanding of 
assembly level language is necessary to understand the developed 
software presented in Appendix F [7]. This microprocessor 
program listing includes the MCS-48 instructions in hex code 
and literal mneumonics and includes full documentation to facili- 
tate interpretation. 

A full performance algorithm was not implemented on the 
Prompt-48 because of its memory storage limitations. The 
original intent was to exercise the software of the complete 


A-7E Takeoff Ground Roll Distance algorithm on the Prompt-48. 
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After the necessary routines were written and stored, only room 
for three coefficients remained (98 coefficients required for 
this algorithm). Since implementation capability was the 
desired result, the computation of a second order polynomial 

was considered sufficient. Although this effort was software 
oriented, the necessary RAM storage for the additional coeffici- 
ents and executive routine could have been easily provided for 

a SBC. 

The software development for algorithm implementation 
Bequired routines for input/output (1/0), executive direction, 
Dinary to binary coded decimal (BCD) and BCD to binary conver- 
mrens, and floating point binary addition and multiplication 
routines. The I/O and executive routines were written by LCDR 
Mreeel. The nonavailability of a number oriented microprocessor 
at the time of this effort required the development of the 
mathematical package described above. The advantages for such 
a capability will be discussed inthe following section. 

In addition to the microprocessor software developed by 
the author and LCDR Siegel, the I/O and display routines would 
require alteration for SBC implementation since a digital dis- 
play and keyboard would replace the Prompt-48. 

Mmeamne 6 Jilustrates the solution method. Figure 7 is 
a flow chart of the program execution sequence. Figures 8 and 9 
Show the Prompt-48 RAM and resident register memory, respec- 


tively. 
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SBC SOLUTION METHOD 


START 


RECEIVE NUMBER AS IT IS ENTERED 






AND PLACE IN BUFFER AS DECIMAL 


DISPLAY BUFFER 


CONVERT NUMBER TO BINARY IN FLOATING 


POINT STORAGE CONVENTION AND STORE 


COMPUTE ALGORITHM 


CONVERT ANSWER BACK TO BINARY CODED 


DECIMAL AND PUT BACK IN BUFFER 


DISPLAY BUFFER 


Figure 6 


Speco lurion Method 
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SBC PROGRAM EXECUTION SEQUENCE 


INPUT INDEPENDENT VARIABLE ''X" 
AND DISPLAY 


CONVERT X TO BINARY NUMBER 
AND STORE 
PLACE By INTO Y-LOCATION 


MULTIPLY (X-LOCATION)*(Y- LOCATION) 


STORE ANSWER IN X-LOCATION 


MOVE By INTO Y-LOCATION 


ADD (X-LOCATION) TO (Y-LOCATION) 


MOVE ANSWER (Byx+B,) INTO X-LOCATION 
[MOVE X_INTO Y-LOCATION | 


MOVE ANSWER (B)x°+B,x+Bg) INTO Y-LOCATION 
CONVERT TO BCD 


Figure 7 


SBC Program Execution Sequence 
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ADDRESS 


000-069 
06A-06F 
Q70-079 
07A-0C6 
0C8-0E2 
OE5-OFF 
100-2EC 


300-3FF 


RANDOM ACCESS MEMORY MAP 


USE 
INPUT AND DISPLAY 
EXECUTIVE ROUTINE SEGMENT 
COEFFICIENT STORAGE 
MAIN EXECUTIVE ROUTINE 
BINARY TO BCD EXECUTIVE ROUTINE 
MISCELLANEOUS SUBROUTINES 
ADDITION AND MULTIPLICATION SUBROUTINES 


BCD TO BINARY EXECUTIVE ROUTINE AND CONVERSION SUBROUTINES 


Figure 8 


Random Access Memory Map 
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ADDRESS USE ADDRESS 
20 LSB 30 
21 X- LOCATION Si 
22 ARITHMETIC REGISTER 32 
23 MSB : 33 
24 EXPONENT 34 
25 |ISB Y-LOCATION | 35 
26 |MSB ARITHMETIC REGISTER! 36 
2) EXPONENT 37 
28 |LSB  BCD-BINARY | 38 
29 IMSB CONVERSION 39 
2A EXPONENT 3A 
2B 3B 
2C 2c 
2D 20: 
2E 3E 
2F 3F 

Jeu ferhats  S) 


RESIDENT REGISTER MAP 


LSB 


DISPLAY HEX 


BUFFER 


MSB 


DECIMAL POINT MASK 


—— eee 


CHARACTER COUNTER 


a me 


LSB 


DISPLAY 


PATTERNS 


MSB 


Resident Register Map 
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The second order polynomial, y = B(0) + BCl)x + B(2)x?, 
was calculated using a mathematical executive routine (alter- 
able for any size polynomial and any number of polynomials). 

The only mathematical operations required were multiplication 
and addition of positive or negative numbers. For speed, binary 
arithmetic was used. For increased storage capability and 
mathematical efficiency, a floating point capability was in- 
eluded. 

The calculation routine proceeded as follows: 


BA "x = (B,x) 


(BAX) a: (Bi) (BA x + BL) 


ss = 2 
(BAX + B) x (Bx + Bx) 


2 : 2 
(Bx “ B,x) + Bo =. (BAX + By 


Although all mathematical operations are performed in 


x + By? 


the 8-bit (2-byte) accumulator register of the 8035 ALU (for 

a SBC, the 8048 PROM), a working accumulator using five regi- 
sters (resident memory registers two through six), was 
established. All numbers in the program (independent variable 
"x" after conversion to binary, coefficients stored in RAM 
070-079 and the 'result') were in one of two binary conventions. 
While in storage, the numbers were in "storage" convention. 

Mie numbers were shifted from "storage" to "working" convention 
only when transferred from the X and Y locations (see resident 
register memory map, Figure 9) to the working accumulator 


(registers two through six). When the desired operation was 


completed, the result was returned to the "storage" convention 
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and moved to the "X" location. Figure 10 displays the "storage 
and working" conventions. 

This software was successfully demonstrated on the 
Prompt-48. The user instructions for the Prompt-48 to repeat 
the demonstration are listed below: 

(1) Ensure the 8035 ALU or 8048 PROM is inserted in the 
"execution" socket of the Prompt-48. 

(2) Enter the program in hex code in the proper storage 
locations as listed in Appendix F. 

On the Prompt-48, press the following keys to clear 
the resident register memory: 
Tel 
“Registers " 


et 


fr 6 
> 


bien! 
tt 8 1f 
Do not press "Program Memory" instead of "registers" 
or the program just entered will be erased. 
(3) To execute the program, press the following keys: 
a 
tf D tt 
"Execute" 
Ge. 
"No Break" 
L1G) tf 
"Execute" 
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BINARY CONVENTIONS 


Each large block depicts 1 byte which includes 4 bits. The compart- 


mented blocks represent 1 bit. 


STORAGE CONVENTION 


' a a oo ie Arn eae aS ae a re 15 WC 
EXPONENT ]—» DECIMAL 
SIGN BIT 


In storage convention, the mantissa is left justified to bit 5. A 
positive number is denoted by 0 in the first bit of the second byte(sign 


bit); a 1 indicates a negative number. 


WORKING CONVENTION 
Me ENE 
2 3 [* 5 «> F F Ite iz 16 12 1% 19| 


(Tet) GEE Gi) Ge 


EXPONENT Ce» DECIMAL 


in working convention, the mantissa is left justified to bit 4. 
The sign bit is stored in FO(X-location number) and F1(Y-location num- 


ber) flags of the program status word. 


Figure 10 


Binary Conventions 
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hMiencdicplaveail blank, awaiting the input of the 

independent variable "x". To enter "x", enter the digit keys 
for numbers (base 10) and "D" for decimal point. "x" will be 
displayed on the digital display as it is entered. To compute 
the algorithm (second order polynomial), press "E". The answer 
Will rapidly appear. To calculate the polynomial with a new 
meee for "x", start at Step 3. 

(4) To prevent the time consuming reloading of the pro- 
gram, it is advisable to store the program on a peripheral 


device (paper tape, disc, etc.). 


2. Single Board Computer using Number Oriented Micropro- 


cessor 

Very recently, the National Semiconductor Corporation 
began production of a chip intended for use in number processing 
meplreations [8]. This chip, the MM57109 MOS/LSI, is capable 
Memeeeescientific calculator functions, test and branch capa- 
bilities, internal number storage, and I/O instructions. Of 
@epecitic calculator functions, only addition, subtraction 
and multiplication would be used. 

A SBC using this chip would need the 8048 PROM for 
coefficient and executive routine storage but would not need 
the space consuming mathematical package of the SBC in the last 
Section. A program counter would still be required but external 
RAM would not. The computation time would be increased over the 
demonstrated SBC (approximate computation time of a HHC), but 
the simplicity of programming would make this proposed SBC very 


Btcractive. 
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wee f-/r TACTICAL COMPUTER 

In February 1978 the author made a trip to the Naval Air 
Facility at China Lake, California. The purpose of this visit 
Meco receive indoctrination on the TC-2/2A tactical computers 
and obtain a programming manual for these devices. The desired 
goal was implementation of selected performance algorithms on 
the laboratory bench computer run by the A-7 Program Office of 
the Naval Weapons Center (NWC). A thorough understanding of 
the computer's capabilities and limitations was provided by Mr. 
Robert Westbrook, a software technician. 

The A-7E computer provides very accurate navigation and 
weapons guidance capability. The TC-2 and TC-2A computers are 
a generation apart, the TC-2A being over two times faster and 
having twice the storage capability of its earlier version. 
BOtn Computers are operational at this time. Specific design 
and programming information is available from the programming 
manual [9]. 

Poe instruction set of the tactical computer provides fixed 
Memmi arithmetic, logical transfer of control (branching), 
address modification and single word input/output instructions 
specifically intended for operations primarily involving arith- 
metic. These features made the implementation of algorithms a 
logical decision. Several factors made this implementation by 
the author impractical. The computer design was quite old, the 
instruction set being very tedious and difficult to interpret. 


Bee computer's inability to function using floating point 
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arithmetic would require a significant software effort in that 
area alone. The time required to become fully familiar with 
the instruction set, write the software, and load and test the 
programs at NWC would have been prohibitive for this investiga- 
ion. 

It is hoped that the programmers at NWC will be able to 
implement those algorithms deemed desireable to achieve an 
onboard capability. Takeoff Airspeed and Maximum Refusal Air- 


Speed are considered ideal for implementation. 


42 





IV. CONCLUSIONS AND RECOMMENDATIONS 


Nine of the A-7E NATOPS Manual performance chart groups 
were reduced to a series of analytical expressions or algo- 
rithms. These algorithms, accurate to five significant figures, 


are as accurate as results obtained by manual manipulation of 











the performance charts. 

Implementation was made on three data processors of different 
programming levels and storage capabilities. These devices and 
degrees of implementation were: 

(1) HP-9830 Desk Computer - complete implementation 

With successful demonstration. 

(2) TI-59 Hand Held Calculator - complete implementation 
With successful demonstration. 

(3) Microprocessor - partial implementation with suc- 
cessful demonstration. 

In view of the success of this investigation, recommendations 
concerning implementation possibilities are listed below: 

(1) Complete reduction of the NATOPS Manual performance 
harts could be accomplished and implemented onto a desk com- 
puter as one large program capable of performance data computation 
ithin seconds. The desk computer would be ideal for mission 
planning on a squadron or air wing level or for Air Operations 
Center use. 

Soe tne programs written for the TI-59 HHC could be con- 


Peetaated onto a CROM and used with a TI-58 HHC for use on a 
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squadron level. As an alternative, the software could be 
rewritten for any HHC of comparable capability. 

(3) Although implementation on a single board computer 
using a number oriented microprocessor is completely feasible, 
because of programming ease and cost consideration, the HHC is 
considered a superior implementation possibility at this time. 

(4) The A-7E tactical computer could easily be pro- 
grammed by software engineers at NWC, China Lake, California, 


to produce an onboard capability. 
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APPENDIX A 


Least Squares Fit Approximation 


References 10 and 11 describe the Least Squares Fit Approxi- 
mation in detail. In general the problem is to represent a set 
of "n'" data points in two-dimensional space 


Viera? 1 born 


Ke 5 
4B ae. 


by a polynomial expression of a curve whose degree is less than 
"nh". Two classes of problems exist: 

(1) Linearly independent - those in which the degree 
(d) of the polynomial is one less than the number of data 
moints 

d = n-1l ip) 

(2) Linearly dependent - those in which the degree (d) 
is less than n-l 

d< n-l (2) 

As an example, a set of four (4) data points randomly spaced 

was chosen. If a third degree polynomial of the form 

ieee BX + CX* + DX? Ce) 
were desired, and the data points x. and Ys were inserted 
(1 = 1 to 4) into four such equations, an exact solution for the 
four unknown coefficients would exist. These four unknowns 
Gould be found from the four equations by numerous conventional 
techniques (Direct substitution, Cramer's rule, etc.). The poly- 
nomial expression generated would be termed a "col-location" 


polynomial because its plot would pass through all data points. 
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It is often advantageous to describe a set of data points 
by a curve that does not pass through each point. This type 
of polynomial would be termed a "regression" equation. For any 
set of data points an infinite number of regression expressions 
exist for any specified degree (een the linearly independent 
case) and the object of the Least Squares Method is to find the 
polynomial coefficients of the chosen degree that best describe 
the data points. In the previous example of four data points, 
assume that, instead of the third degree form chosen, a second 
degree equation were selected of the form 

aa et Be OX (4) 

With four data points, the polynomial 1s overspecified and thus 
linearly dependent. For this case an infinite number of solu- 
tions exist for the coefficients a, b and c. If an error term 
($) were defined for any given X,Y pair as 


ae fe BX 4 Cx. | C5) 


oa] ale 
a total squared error term (E) could then be defined by squaring 
and summing the terms attained: 
N 
E= 2 6,’ (6) 
i=l 
If E were them minimized for any given degree chosen, the best 
Least Squares Fit would have been achieved. 
If the values for 6 from Equation 5 were inserted in Equation 
6 and the partial derivative of E were taken with respect to the 
coefficient A, an equation would be generated that when set equal 


to zero (0) would define a minimum value of E for a given value 


of A. If the same operation were performed with respect to the 
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coefficients B and C then three equations would be generated 
With three unknowns (A, B and C). The solution of these simul- 
taneous equations would produce the coefficients A, B and C, 
that would minimize the value of E and hence would produce a 
Least Squares Fit approximation to a set of linearly dependent 
equations. 

A numerical procedure has been developed to accomplish this 
task. An example of this procedure has been included in the 


following paragraphs [10, 11]. 


Least Squares Fit Method Example 

Given the following set of data: 

x 0 Hi 2 ~ 7 
| BCX) =Y 0 ils 3 12 20 
Meet a Curve of the form 

mG = Y - AU BX + CX? 

STEP 1: Substitute all pairs of data into the form equation 
Mepbelding the fact that the coefficients (A, B and C) must 


Satisfy all the following: 


i= A+ B(O) + COO)? 
dee + B(1) + C(1)* 
oe A+ B(2) + C(2)? 
fe A BCH) + CCH)? 
moe= A + BC7) + C(7)? 


Now multiply each expression by its coefficient of A in that 


expression and add all equation yielding 
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coe] Jf, LY Peres OC 
Now multiply each expression by its coefficient of B in that 
expression and add all the equations yielding 


He=- OCA) + OCB) + OCC) 


is 4 f+ IB + IC 

6 = 2A + 4(B) + 8(C) 
Hs = 4A +16(B) +64(C) 
140 = 7A +44U(B) +343(C) 


moo = LYA +70(B) +416(C) 
Now multiply each expression by its coefficient C in that 
expression and add all the expressions yielding 
wee OCA) + OCB) + OCC) 
ime cA) + 1(B) + 1¢C) 


a2 


heen 8B) +16(CC) 
me2>=16(A) +64(B) +256(C) 
meu -49(A) +343(0B)+2401(C) 
meses = /OA + 4HI6B + 2674HC 
Now solve the following three previously generated equations for 
Mme coefficients A, B and C yielding 
fee = SA + 1YUB + 70C 
195 =14A + 70B +416C 
ies =/0A +tH16B +2674C 
Pees, Fo = 2.6, C = .065 
and 


fae SO +t 256 X% + 2.065X- 
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The following plot and chart depict the original data and 


the data obtained from the equation for the fitted curve: 


Original Fitted Curve 


Polynomial 
Xx nC he 
0 0 esl. 
l i 1.67 
2 3 iene 
u ‘2 10.46 
; AMO 20.41 
Q.E.D 


Bhs) 
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APPENDIX B 


NATOPS Manual Performance Charts 


These charts from which the performance algorithms were 


developed are listed below in order: 


Figure 


Bl 
B2 
B3 
Bu 
B5 
B6 
B7 
B8 


Bg 


B10 
Bid 
ba 2 


bi 3 


Crmrce 
Cruise 
Gminise 


Crisco 


Title 


Performance, 
Performance, 
Performance, 


Performance, 


Takeoff Factor 


Phase I 


Piiese sl 


Phase III 


Phase IV 


Takeoff Ground Roll Distance 


Adjusted Takeoff Ground Roll Distance 


Maximum Range Cruise at Constant Altitude (Time, Speed) 


Maximum Range Cruise at Constant Altitude (Fuel 
Required) 


Military Power Climb Schedule 


Takeoff Speed 


Maximum Refusal Speed 


Cruise Ceiling and Optimum Endurance Altitude 
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Cruise Performance, Phase IV 
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Figure B6 


Takeoff Ground Roll Distance 
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Figure B7 


Adjusted Takeoff Ground Roll Distance 
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Figure B9 


Maximum Range Cruise at Constant Altitude (Fuel Required) 
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Figure B10 


Military Power Climb Schedule 
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TAKEOFF SPEED (A-7E) 


ENGINE: TF41-A.2 
FUEL GRADE: JP.5 
FUEL DENSITY: 6.8 LB/GAL 





MILITARY RATED THRUST 
LANDING CONFIGURATION 
LEADING EOGE FLAPS OTN 


CONDITIONS: 


DATA BASIS: FLIGHT TEST 
OATE: NOVEMBER 1971 


MODEL: A-7E 


Increase speed 1/2 knot per 1% 


Decrease speed 1/2 knot per 1% aft 


Oate basis is 26% MAC, 
forward CG shitt. 


CG shift 
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CRUISE CEILING AND OPTIMUM ENDURANCE ALTITUDE (A-7E) 





MODEL: A-7E ENGINE: TF41-A-2 
DATA BASIS: FLIGHT TEST FUEL GRADE: JP-5 
DATE: NOVEMBER 1971 FUEL OENSITY: 6.8 LB/GAL 
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Cruise Ceiling and Optimum Endurance Altitude 
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APPENDIA C 


Generated Algorithms 
BOW LEVER CRUISE PROGRAM 


Pnase [ 


[Ml = -92.512 + 236.8966 


| 


‘Transfer scale Versus Drag Count 

Me 2.3287 -.26316D + .0073327D2 -(7.513E-5)D® + (3.5396E-7)D* 
-(7.78E-10)D° + (6.462E-13)D° 

Memisss + 1.0956D -.030653D2 + (3.1912E-4)D? -(1.5276E-6)D" 
+(3.408E-9)D° -(2.8692E-12)D° 

m2 = 10.284 -1.0719D + .031094D? -(3.2878E-4)D°* + (1.595E-6)D* 
-(3.6009E-93)D° + (3.0634E-12)D° 


mme= AQ + (Al1)(M1) + (A2)(M1)? 


Transfer Scale Versus Guidelines 

fee= 22.819 -31.7341I + 41.3317 -5.0953I° 
Me =-154.98 + 217.512 -261.7317 + 35.905I1° 
B2 = 405.08 -525.56I + 607.49I* -88.7371° 
Me = -445.62 + 542.981 -611.5517 + 92.8941" 
Me 104.78 -204.42T + 225.891% -35.1891° 


Me= BO + (B1)(M1L) + (B2)(M1)2 + (B3)(M1L)° + (B4)(M1)* 


Mmaase IL 


Mees +t 2(6(4.3732E-3) + .027743D]M? 
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Phase IIil 


BO 


ome 


B2 


B3 


ee 


R1 


2 


N1 


N2 


6.6253 -1.989R + 3.0252R* -1.0761R° + .17675R” -.013095R° 
+ (3.526E-4)R° 

MiiercOl2 —248.9317R * 91.66355R* -15.55218R° + 1.224432R" 
-~.0395333R° + (2.896385E-4)R® 

S22 123 + 1231.24R -487.4233R* + 91.6522R° -8.662962R* 

+ .3953974R° -~.006905535R° 

1680.142 -1950.139R + 788.8513R* -152.5733R* + 15.03819R°* 
mee? 274139R° + .013707R° 

R 

2 (Integer (R/2)) 

Rl + 2 

ple (B1)(R1) + (B2)(R1)? + (B3)(R1)° 


Bore (B1)(R2) + (B2)(R2)2 + (B3)(R2)? 


Using Linear Interpolation 


een! + ((N2-N1)(R3-R1)/2] 
P = 4,974U6N + (7.9043E-6)N? 
Biase LV 
Memeo. 4375 + .010426T -(6.8925E-6)T* + (€4.9127E-7)T?* IM 
m= .LC(N4Y)P 
TAKEOFF DISTANCE PROGRAM 
Meme a.086 —-.00017113A -(2.0655E-7)A2 + (3.6861E-11)A’ 


mye? 4156E-15)A° 





| Bl 


B2 


IBS 


ie 


ia 


GO 


Gl 


ist 


If 


LO 


L1 


ot 


men g45635 —-(7.8931E-6)A + (3.7545E-9)A* -(9.7088E-13)A°? 
mec Gg .997E-17)A‘* 


mol? =C8e2558E-7)A + C4.0739E-10)A* -(8.548E-14)A° 


+ (5.4964E-18)A" 
= ~(1.9097E-5) + (1.3671E-8)A -(9.4694E-12)A* + (2.0434E-15)A° 


mie 617E-19)A° 


BO + B1(B) + B2(B)? + B3(B)? 
Boupile datum on, 

Smee? 73 + .56598C 

double datum off, 


meoyi7s + .65876C 


= -(4.8896E+5) + (8.4974E+1)G -(5.7856E-3)G* + (1.9373E-7)G° 
-(3.1744E-12)G* + (2.0446E-17)G° 

= (5.8621E+4) -(1.0146E+1)G + (6.8807E-4)G* -(2.292E-8)G° 
+ (3.7387E-13)G* -(2.3964E-18)6° 

=)G0 + G1(E) 

relative humidity < 40%, K = H 

not, K = 4¥{{(I-40)/1000]+1} 


meeou + .89509K + (2.3306E-5)K* -(1.6254E-9)K? 


tl 


fees. 3728E-14)K* 


Memes) —(il.08566-2)K + (2.1754E-7)K* -(2.5327E-11)K’? 


mel, 197E-15)K' 


Mie? = -1) —-(?.16666-6)K + (3.4274E-9)K* -(2.7817E-13)K° 


fmm, 00 77E-18)K* 
melo + L1(L) + L2(L)? 


Pandas calm, M = K 
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XO 

Xl 
2 
—Q0 
Ql 


30 
ol 


#4 


BO 


(4.5704E+1) + .93429M + (2.2265E-5)M* -(2.338E-9)M? 
yer =Lh Mm 

7.9472 + .O14914M + (9.0708E-6)M* -(7.1235E-10)M? 

+ (3.0684E-14)M" 

5.3616 -.0085136M + (3.5914E-6)M* -(4.5932E-10)M? 

+ (1.9889E-14)M’ 

ee; X(N) + X2(N)? 

@604%.2 -2.1694X + .0010915X* -(1.1119E-7)X* + (3.662E-12)x’ 
sumo 73 + .22601X -(7.5225E-5)X* + (7.7018E-9)xX? 
-(2.5437E-13) xX" 

moe so =. 0040102 + (1.2832E-6)X* -(1.3234E-10) xX? 

+ (4,3908E-15)xX"* 

QO + Q1(P) + Q2(P)? 

mover 79) + 1.58010 -(2.0254E-4)Q? + (2.4111E-8)0° 
-(8.6737E-13)Q" 

16.196 -.024333Q + (9.3484E-6)Q? -(1.2594E-9)Q? 

+ (4.7522E-14)Q* 

meeerss + (2.359E-4)Q -(1.037E-7)07 + (1.6016E-11)0° 
=(6.3195E-16)Q" 


eie+ SLC(R) + $2(R)? 


Pe eve ReaniGem eh Ulor) PIME AND SPEED 
AT CONSTANT ALTITUDE PROGRAM 


Pmt (5.0794E-3)H -(1.3968F-3)H* + (8.254E-5)H? 


men. 2698E-6)H* 
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Bl 


BO 


Bl 


B2 


B3 


M1 


‘al 


BO 


Bi 


B2 


.05 + .0015159H + (1.123E-4)H* -(3.4921E-6)H? 

+ (7.9365E-8)H”* 

BO + B1(G) 

547803 + .0013417D + (6.2287E-6)D* ~(1.6261E-8)D* 

+ (1.6438E-11)D* 

.08217 + (4.1209E-4)D -(4.5577E=-6)D? + (1.6777E-8)D° 
-~(2.001E-11)D’ 

(4.2143E-4) -(9.4397E-5)D + (1.2646E-6)D* -(4.8537E-9)D° 
fms. 7222E=-12)D' 

-(6.6767E-4) + (8.4671E-6)D -(1.0501E-7)D? + (3.6382E-10)D° 
sis). 7828E-13)D’ 

BO + B1(N) + B2(N)? + B3(N)? 

M -{(60-T)(2)(M)/1200] 

(710)(M1-.14) + 100 -E 


D1/V 


PUPoe eOUlrREDOTOR MAXIMUM RANGE CRUISE 
AT CONSTANT ALTITUDE PROGRAM 


4U.54 =~. 16444UA + .0033932A* -(1.0283E-4)A* + (1.926E-6)A* 
=(1.3757E-8)A° 

(3.22E-9) -(3.6664E-3)A + (8.9338E-4)A* -(5.5939E-5)A° 

+ (1.4593E-6)A* -(1.3281E-8)A° 

(6E-4) + (1.1203E-4)A -(2.3358E-5)A* + (1.4536E-6)A° 
Bs. 7144E-8)A* + (3.3334E-10)A° 


Paet B1(G) + B2(G)? 


68 





BO 


Bl 


Ul 


V0 


V1 


V2 


V3 


BO 


Bl 


-(2.5399E-3)D + (9.7299E-5)D* -(2.3516E-7)D? 

fc. 42515-10)D" 

2 + (4.2388E-3)D + (1.2326E-5)D*? -(1.0298E-7)D° 
fo. 7277E-10)D" 

BO + B1(N) 

ey 


@.)(1)/60 
Pe oni kewiee, CLIMB ATRSPEED SCHEDULE 


Mmiemeso —.79075D + .0011382D* -(4.1018E-7)D° 
meme et 6 24E-3)D + (7.6582F-5)D* -(1.1344E-6)D? 
meme 71295E-9)D* -(2.3035E-11)D° + (€3.6588E-14)D° 


Pre? 3062E-17)D’ 
TAKEOFF AIRSPEED PROGRAM 


54.023 + (3.4787E-3)G -(1.9475E-8)G? 

ae ((26-P)/2] 

seeo7. 1 + 61.604U -.70348U2 + .0035661U°* -(6.6578E-6)U" 
Meme 2s =2.4517U + .028779U* -(1.4753E-4)U* + (2.7872E-7)U' 
=.72239 + .023415U -(2.798E-4)U* + (1.4596E-6)U? 
~(2.807E-9)U"* 


VO + VICR) + V2(R)? 
MAXIMUM REFUSAL SPEED PROGRAM 
mmo) + 6.761G —-.35159G* + .0080545G° -(6.7769E-5)G° 


26.312 -3.8382G + .20326G? ~.047022G* + (3.994E-5)G" 
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iB 2 


B3 


BO 


Bl 


B2 


=4.9639 + .72723G ~.038721G* + (8.985E-4)G? -(7.638E-6)G" 
30288 -.044855G -.0023921G* -(5.5549E-5)G? 

+ (4.7217E-7)G* 

BO + BL(E) +B2(E)* + B3(E)? 

-~11.412 + 62.185L -9.0037L2 + .64921L° -.017455L" 

Beye) -§ .2012L + .70377L* -.058693L° + .OO1L7461L* 


BO + BLCR) 
Op ium ENDURANCE ALTITUDE PROGRAM 


55.333 + .073076D -(9.7836E-4)D? + (3.5015E-6)D? 
-(3.9782E-9)D* 

~1.1 -(8.0597E-3)D + (8.0097E-5)D* -(2.8836E-7)D° 
es. 3032E-10)D" 

(6.6667E-3) + (1.2541E-4)D -(1.4039E-6)D?* 


pom BL(G) + B2(G)? 
CRUISE CEILING PROGRAM 


85.118 -.29117D + .0030434D? -(1.2851E-5)D* + (1.6621E-8)D* 
Bee7s77 + .025635D -(3.3063E-4)D* + (1.4162E-6)D° 
-(1.8343E-9)D* 

§063327 -(8.5289E-4)D + (1.0814E-5)D* -(4.6514E-8)D° 

+ (6.0606E-11)D* 

=(6.0468E-4) + (9.0826E-6)D -(1.143E-7)D* + (4.9304E-10)D° 
-(6.4567E-13)D* 


BO + (BL)G + (B2)G? + (B3)G° 





- APPENDIX D 
HP-$830 Programs and Lists of Variables 
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H2=H 
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RET COMPLETED CALCULATION OF INTERMEDIATE # BY LIHERF INTERFOLATION 
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PRIHT “GROSS HT="G+ Loe 

PRIHT “Tss"S°0C="D°M="h 

FRIHT “TEMF="T 

PRIHT “REF #="F3 

PRIHT “H="H 

PRIHT “LBFIUEL- Hits" F 

PRINT “FUEL FLOW= "F 

PRIHT 

BOTH te 


eee eo tee let th sae LP e-5. 9534103 
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List of Variables for Program 1 


Variable Henin. Lon 
G Gross weight (1bs..) 
D Drag count 
als Temperature (°C) 
M Mach number 
M1 Result of lower graph, Figure Bl 
18 Guidelines, numbered top to bottom consecutively 
> Transfer Scale calculated as function of I 
1 iMaictem scale Calculated as function of D 
a2 Transfer Scale calculated for upper guideline 
38 Transfer Scale calculated for lower guideline 
Tt Relative Transfer Scale location between guidelines 
fR3 Reference number 
Ri Even reference number below actual reference number 
i 2 Even rererence number above actual reference number 
J Integer counter 
N Result of lower graph, Figure B3 
Nl Result of lower graph, Figure B3 for Rl 
N2 Result of lower graph, Figure B3 for R2 
N4 Result of lower graph, Figure B4 
a o Coefficients 
. Coefficients 
iP Pounds of fuel per nautical mile 
rE Fuel flow 
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Program 2 


PEM THIS PROGRAM CALCULATES THE TAREVUFF DTSTAHLE FEWUIFED FOR AW A-7E 
FEM IT IS DEFENUENHT OH 3 VARIABLES -- 
FEM GROSS WEIGHT-«FHUY ALTITUDE: TENE: DRAG COUHT RELATIVE HUMBDITY  sWIHDS 
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feet =THREGQFF ROLL DIST="S 

LOTO 3 

EHD 


74 





Variable 
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D 
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List of Variables for Program 2 


Derinitton 

Runway Altitude (feet) 

Temperature (°C) 

Detiremcdatume status (1 indicates "with") 
Gross weight (1bs.) 

Relative humidity (%) 

Headwind (kts.) 

Runway slope (%) 

Comter Of Gravity (3) 

Flap position (degrees) 

Result of upper graph, Figure B5 

Takeoff factor 

Mnedyusted ¢round roll distance, Figure B6 
Adjustment factor due to relative humidity 


Ground roll distance (GRD) adjusted for relative 
neanGl Jae y- 


GRD adjusted for wind 
GRD adjusted for runway slope 
GRD adjusted for the center of gravity location 


True GRD (also adjusted for flap position) 


Coefficients 
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Coefficients 
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Program 3 


FEM THIS PROGRAM CALCULATES THE A-TE MASINI FRHGE Al 
Sere TIME OF FLIGHT AHD 15 DEPENDENT Ut « YARREIRELES - 
PEM GROSS WEITGHT« ALTITUDE: OFAG COU). TENFERRTUREs WIND: 
Seey))6 CULHPUT 6G ALT: D0 TEMP eo HOMO. ODISTRAHCE” 
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List of Variables for Program 3 


Variable Definition 
G Gross weight (l1bs.) 
H farex rude (ft. ) 
D Drag count 
iT Temperature (°C) 
L Headwind (kts.) 
D1 Distance to fly 
N Rest@peon first chart, Figure B8 
M Cruise Mach number (adjusted and unadjusted for T) 
V Ground speed (kts.) 
al ome of flaght 
A0 ,BO Coefficients 
Pa, BL @eetficients 
B3 Coefficient 
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REM 
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FR INT 
PRINT 
LHEUT 
FR IIHT 
FR INT 
FR IHT 
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FR IHT 


ALTITUDE FOr AH A- 
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Program 4 


THIS PROGRAM CALCULATES FUEL FPEGUIRED FOR MAM FAHGE AT CONSTANT 
rE AHD IS DEPENDENT OH OS VARIABLES 

HE TGA eee PIPE DRAG COUNT: TRUE AHIFSFEEDs AnD 

“EMTER GROSS Ws ALT+ DRAG CT TAS: TIMES MINUTES 2" 


ole ani) ec, 
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List of Variables for Program 4 


Variable Definition 
G Gross weight (lbs.) 
A Altitude (ft.) 
D Drag count 
V True airspeed (kts.) 
T iemeson tlvent (minutes) 
N Resiiupeonearst chart, Figure B9 
L Poumads of tUel per nautical mile 
iz Fuel flow 
R Fuel required 
BO ,A0 Coefficients 
B1l,Al Coefficients 
B2 Coefficient 
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Program § ‘ 


1 FEM THIS PROGRAM CALCULATES THE CLINE AIRSPEED OF AH A-Fe 
> REM ITHDICATED AIRSPEED BELOW 2o.008°" » 
REM GACH HUMBER ABOVE 26s? 3 


MW poe 
S PRIHT “ELINE AIBPSPEED SCHEDULE” 
Sereind! “DRAG Lt CLINE AIRSPEED EEE Hae 


Poe Oe a 


mereIHT “ Sse EE Shes 9” 
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42 M=Mel a 
44 M=THT sts 
$4 A eco 
MeerelHy Des; it 
be DS=D+1 sy 
Seer Ossi THEN 24 
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List of Variables for Program 5 


Dement Om 
Drag count 


Mach number 


Calibrated airspeed (kts.) 
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Program 6 


Sees FRUGEAM LALLULATES THE TAKEOFF AIRSPEED OF AN A-TE 
REM UWHDERF VARYIHG GROSS WEIGHTS: FLAP POSITIONS 

See AHI CEHTER OF GRAVITY LULATIONS 

Rec 

Bee" 


PRIHT "ELIP BROS WETGHT="G 

FR IHT 

PRIMT 

PRIHT “FLAPS CG TAKEOFF AIRSPEED” 
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Variable 


UL 


V4 


7 Vl, 
V2,V3 


List of Variables for Program 6 


Deere on 

Flap position (degrees) 

Center of gravity (%) 

Gross weight (l1bs.) 

Unadjusted takeoff airspeed 

Takeoff airspeed adjusted for center of gravity 


Actual takeoff airspeed (adjusted for flap position) 


SCoOCrrICLents 
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Program 7 


Bement FRUVNGRAM CALCULATES THE NAS TMUM REFUSAL SFEEB 
REM FOF AH A-rE USI AHTI-SEID 
eee) Uk 6 PEFPEHDEHT OH S VARIABLES -- 
FEM GROSS HETGHTs TEMPS FHM Y LCEHGTH FHMY ALTITUDE: AHO DOUGBLE DATUM 
PETHT “THPUT ALT: TEMPS RMN CTHs GH,» DOUELE DATUM" 

IHFUT As GslsGe0 

B=5- 16h 
L=L“1tek iat. 
Ser) 3=60AL TI TU0E="A 
PRIHWT “TEMNF="E 
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FR THT “TRIN WTS" ele 
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ge SSL is SUI SEL ti. Tes Pelt o-b. asses seL tse. bol Pde l elt 
MN M=BO+61 +k 

M=sTHT om 

PRIHT . 
Seer ttl * TAREGQFF FACTOR ="E 
eee TIT 
eee) MH FEFUSHL SPEEB = “st 
1 FRIWT 
MH GOTO F 
HM EMD 
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Variable 


A 


8 i 2 a Se Si ae 


Be, Bl, 
Be, 83 


hiSewon Vartvaples for pare ars ual oi 


Ise sig aL Te abengl 


Runway Altitude (ft.) 

Temperature (°C) 

Runway length (ft.) 

Gross weight (lbs.) 

Double datum status (1 indicates "with") 
Result of upper chart, Figure BS 

Pakeorft factor 

hesult Of first chart, Figure B12 


Maximum refusal speed (kts.) 


Coefficients 
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Program 8 


FEM THIS PROGRAM CALCULATES THE OPTIMUM EMHDURAHITE ALTITUDE 
REM OF AH A-°E AT VARYING GROSS WEIGHTS AnD DRAG COUNTS 
err EC: ] 

Gb=13 

[=6 

PRIHT “OPTIMUM EMHDURAHCE ALT 

PRINT SGeitss Wy DRAG CT WET EHO ALT” 
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Mm H=BC3I+60 1 J+G+El 2 I-Gr2 

Seeo=(H) Hei weus 
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em IF BL 318 THEN 28 
¢ [=8 

3 IF G.45 THEN Se 
M EMD 


86 





List of Variables for Program 8 


Variable Detinxcion 
G Gross weight (lbs. times 1000) 
D Drag count 
H Optimum endurance altitude (ft.) 
i Optimum endurance altitude ( integer format) 
x Gross weight (lbs.) 


paeB2,B3 Seerfricilents 
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Program 9 
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EusemOuevadriaples for Program 9 


Variable Befination 
G Gross weight (lbs. times 1000) 
D Drag count 
H Gase cealang (ft.) 
a Cruise ceiling (integer format) 
X Gross weight (l1bs.) 

si Coefficients 
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TI-59 Programs and User Information 


USER INFORMATION FOR PROGRAM 1 
Program: Low Level Cruise Performance 


Number of Steps: 1386 


Computation Time: 90-110 seconds 
STEP ENTER PRESS KEY 
1 gross weight (l1bs.) A 
2 seas count C 
c mach number D 
4 temperature (°C) iS 
5 "--- R/S 
6 read in cards 3 & 4 - 
7 drag count C 
8 mach number D 
9 temperature (°C) E 
0 ee nie 
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Dis Picea 
gross weight/1000 
diese coum 
mach number 
Transfer Scale 
Unusable number 
Transfer Scale 
mach number 
lb.fuel/nautical mile 


fuel flow 
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USER INFORMATION FOR PROGRAM 2 
Program: 
Number of Steps: 


Computation Time: 
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gross weight (lbs.) 
pressure altitude (ft.) 
temperature (°C) 
headwind (kts.) 
relative humidity (%) 


read in cards 3 & 4 
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center of gravity (%) 


runway slope (3%) 


48-50 seconds 


Takeoff Ground Roll Distance 
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USER INFORMATION FOR PROGRAM 3 


Program: Maximum Range Time and Speed at Constant Altitude 


Number of Steps: 547 


Computation Time: 20 seconds 

STEP ENTER PRESS KEY 
BL erPoss weight (libs.) A 

2 peltatude (ft.) B 

3 drag count C 

u temperature (°C) E 

5 distance (nautical miles) 2nd,B 

6 headwind (kts.) Zool at 

7 --- RS 

8 --- sos 


aecl2 


DISPLAY 
gross weight/1000 
altitude/1000 
dra treo Unme 
temperature 
distance 
mach number 
ground speed (kts.) 


time required (minutes) 
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USER INFORMATION FOR PROGRAM 4 
Program: Maximum Range Fuel Required at Constant Altitude 
Number of Steps: 439 


Computation Time: 16 seconds 





STEP AN ee PARC Sys) ag DISPLAY 

1 gross weight (lbs.) A gross welght/1000 

2 pubtatude (ft. ) B altitude/1000 

5 aeae Count G drag count 

n true air speed (kts.) PargheliaD) true airspeed 

5 time (minutes) /Laala ie (lb.fuel required/ 
nautical miles) 

6 --- RS fuel flow required 

7 --- Ryo fuel required 
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Moon INFORMATION FOR PROGRAM 53 


Program: 


Number of Steps: - 


Computation Time: 


BL LP 
aM 


2 


ENTER 
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USER INFORMATION FOR PROGRAM 6 
Program: Takeoff Airspeed 


Number of Steps: 2605 


computation Time: 93 seconds 

STEP ENTER PRESS KEY DISPLAY 

lL gross weight (lbs.) A gross weight 

y flap position (degrees) 2nd,A flap position 

3 center of gravity (%) Pia. B takeoff airspeed (kts) 
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USER INFORMATION FOR PROGRAM 7 
Program: Maximum Refusal Speed 


Number of Steps: 716 


Computation Time: 26 seconds 
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1 gross weight (1bs.) A 

2 runway altitude (ft.) B 

3 temperature (°C) i 

4 runway length (ft.) Za 
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USER INFORMATION FOR PROGRAM 8 
Program: Optimum Endurance Altitude 


Mumber of Steps: 229 


Computation Time: 8 seconds 
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USER INFORMATION FOR PROGRAM 9 
Program: Cruise Ceiling 
Number of Steps: 299 


Computation Time: 11 seconds 


BLEEP ENTER eS Saapd DISPLAY 
i gross weight (lbs.) A gross weight/1000 
2 drag count G cruise celling 
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Neve 1a; DAL.G 13 


Microprocessor Program 


This computer printout was reproduced on an IBM-360/67 
Computer. This listing was entered by hand on an INTEL Corpora- 
tion MDS-800 System, stored on a diskette, and sent to the 


IBM-360/67 via a high speed line. 
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